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The  structure  of  echoes  observed  in  echo-location  systems  such  as  radar  and  sonar  can  exhibit  significant 
variability  caused  by  environmental  effects  and  fluctuations  in  the  target?s  aspect.  A  technique  for 
countering  these  effects  and  thereby  improving  the  detection  and  classification  of  the  target  is  investigated. 
The  technique  relies  on  being  able  to  characterizes  the  target?s  signal  space  using  a  low  dimensional  signal 
basis.  For  this  purpose  the  Karhunen-Loeve  expansion  is  used  to  constmct  the  target?s  signal  space  and 
in  conjunction  with  the  likelihood  ratio  test  form  the  appropriate  test  statistic,  which  is  the  basis  of  the 
detector/classifier  stmcture  for  distinguishing  target  returns  from  nontarget  ones.  This  structure  is 
sometimes  referred  to  as  the  separable  kernel  receiver.  With  this  stmcture  in  place  the  performance  of  the 
detector/classifier  can  be  evaluated  as  a  function  of  clutter  type  and  target  mismatch.  Performance 
comparisons  are  also  made  with  a  peak  detector  and  a  matched  filter  processor  detecting  Rayleigh  and 
Rician  targets  in  Rayleigh  clutter.  The  results  are  presented  in  the  usual  form  of  receiver  operating 
characteristic  (ROC)  curves. 
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DETECTION  AND  CLASSIFICATION  OF  FLUCTUATING  TARGETS 


1  INTRODUCTION 

The  detection  and  classification  of  targets  using  echo-location  systems,  such  as  sonar  or 
radar,  can  be  difficult  depending  on  a  number  of  factors.  Two  of  the  most  important  are  the 
environment  in  which  the  target  operates  and  the  target  itself.  For  example,  the  environment 
can  produce  adverse  propagation  and  clutter1  to  obscure  the  target,  while  the  target  itself  can 
present  many  different  physical,  and  hence  acoustic  or  electromagnetic,  profiles  to  the  echo- 
location  sensor,  making  its  echo  potentially  difficult  to  detect  and  identify  on  a  consistent 
basis. 

Gaumond  in  a  recent  paper  [1]  presented  a  technique  for  the  detection  and  classification 
of  a  hemispherically  end-capped  cylindrical  shell  with  periodic  internal  ribs.  The  essential 
idea  was  to  show  how  deterministic  signals  from  a  model  of  target  scattering  can  be  used 
to  develop  a  detection/classification  scheme  for  a  target  that  has  an  aspect  dependent  echo 
response.  In  general,  the  aspect  angle  of  the  target  was  unknown  and  assumed  to  be  randomly 
distributed.  Therefore  it  was  appropriate  to  model  the  target  echo  as  a  stochastic  signal.  This 
signal  was  then  decomposed  into  a  set  of  eigenfunctions  using  the  Karhunen-Loeve  expansion 
(KLE).  It  was  found  that  the  target  echo  was  well  characterised  by  a  small  number  of  KLE 
terms.  These  eigenfunctions  and  their  eigenvalues  were  then  used  to  form  a  set  of  filters  to 
detect  and  classify  the  target.  For  this  particular  case,  the  stochastic  target  response  was 
based  upon  physical  components  of  the  echo,  e.g.,  the  Bloch  and  Chalice  waves  that  are 
readily  evident  below  1  kHz  [2], 

The  basic  idea  of  this  physics-based  approach  appears  to  be  sufficiently  general  to  allow 
for  adaptation  to  higher  frequencies  even  though  it  is  no  longer  obvious  what  are  the  relevant 
physical  features  that  can  be  used  for  detection  and  classification.  The  purpose  of  this  paper 
is  to  extend  the  KLE  technique  to  high  frequencies  where  physical  features  may  not  be 
readily  apparent  and  the  target  is  embedded  in  non-Rayleigh  clutter.  At  high  frequencies 
the  impulse  response  of  the  target  echo  appears  as  a  superposition  of  scattering  highlights 
which  changes  with  its  aspect.  It  is  hoped  that  by  projecting  the  target  echo  response  (as  a 
function  of  aspect  angle)  onto  a  signal  space  representation  via  the  KLE,  a  useful  set  of  “signal 
templates”  of  low  dimensionality  can  result.  In  effect,  the  KLE  acts  as  a  feature  extractor.  If 
this  is  the  case,  then  it  would  be  of  interest  to  see  how  the  detector /classifier  would  perform 
in  a  highly  cluttered  environment.  A  non- Rayleigh  representation  of  the  clutter  is  chosen 
because  of  the  higher  probability  of  getting  large  amplitude  events  mimicking  false  targets. 

Manuscript  approved  December  ??,  2000. 

Hn  this  report,  “clutter”  refers  to  any  signal  that  is  composed  of  a  continuous  echo  from  the  environment. 
This  term  is  commonly  used  in  radar  parlance.  In  sonar  parlance,  the  equivalent  term  is  “reverberation.” 
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The  outline  of  the  rest  of  the  paper  is  as  follows.  In  Section  2  the  derivation  of  the 
test  statistic,  and  therefore  the  detector/classifier  structure,  is  given.  In  Section  3  the  de¬ 
tector/classifier  is  applied  to  a  test  example.  The  target,  signal  and  background  (clutter) 
models  are  specified.  This  is  then  followed  by  a  performance  evaluation  in  Section  4  in 
which  receiver  operating  characteristic  (ROC)  curves  are  presented,  showing  comparisons 
with  the  cases  of  Rayleigh  and  Rician  targets  in  Rayleigh  clutter,  as  well  as  a  comparison 
with  the  peak  detector.  Furthermore,  the  detector/classifier’s  performance  as  a  function  of 
clutter  type  (from  Rayleigh  to  non- Rayleigh)  and  target  mismatch  is  also  given.  Finally, 
some  conclusions  and  comments  are  given  in  Section  5. 


2  DETECTOR/CLASSIFIER  STRUCTURE 


This  section  reviews  the  derivation  of  the  detector /classifier  structure  and  explains  the 
motivation  and  assumptions  behind  the  choice  of  the  structure.  The  derivation  is  based  on 
the  problem  of  detecting  a  coloured  Gaussian  target  in  white  Gaussian  clutter. 

A  standard  and  well  established  way  of  determining  the  structure  of  a  detector  is  to 
derive  the  test  statistic  from  the  likelihood  ratio  of  the  target-plus-clutter  and  clutter-only 
probability  density  functions  (PDFs)  [3].  Ideally,  an  optimal  detector  can  be  derived  if  the 
target  and  clutter  PDFs  are  known.  In  practice,  this  is  generally  not  the  case  as  the  char¬ 
acteristics  of  both  the  target  and  clutter  can  change  spatially  and  temporally  for  scenarios 
of  operational  interest.  In  some  other  cases  the  target  and  background  information  may  be 
incomplete  or  simply  not  known  at  all.  So  instead  of  seeking  an  optimal  structure  it  is  as¬ 
sumed  that,  for  the  purpose  of  deriving  a  useful  test  statistic,  the  clutter  is  a  zero  mean  white 
Gaussian  random  process  and  the  target  is  represented  by  a  zero  mean  coloured  Gaussian 
process.  The  choice  of  statistics  for  the  clutter  is  motivated  by  the  central  limit  theorem  and 
is  thought  to  represent  the  limiting  case  of  many  situations.  Similar  considerations  follow 
for  the  target  but  with  less  restrictions  on  the  form  of  its  correlation  function  (i.e.,  sample 
returns  from  the  target  response  are  in  general  statistically  dependent).  Clearly,  the  per¬ 
formance  of  the  detector  will  be  less  than  optimal  once  these  assumptions  are  violated,  but 
just  like  the  matched  filter,  the  simplicity  of  the  resultant  structure  may  have  the  virtues  of 
robustness  and  utility  for  a  wide  range  of  situations.  The  usefulness  of  such  a  detector  can 
then  be  tested  by  subjecting  it  to  a  number  of  non-ideal  situations  of  interest. 

By  using  the  above  assumptions  the  detector  structure  can  be  derived  as  follows.  To 
begin,  consider  the  binary  detection  problem  with  the  null  hypothesis  ( H0 )  that  no  signal  is 
present  and  the  alternative  hypothesis  (Hi)  that  the  target  is  present.  Symbolically: 

H0  :  r(t)  =  n(t) 

Hi  :  r(t)  =  s(t )  +  n(t)  (1) 

where  r(t )  is  the  received  waveform,  s(t)  is  a  realisation  of  a  zero  mean  coloured  Gaussian 
process  representing  the  target,  and  n(t)  is  a  sample  waveform  of  a  zero  mean  white  Gaussian 
process  representing  the  clutter.  Note  that  all  waveforms  are  assumed  to  be  complex  valued 
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for  generality.  Since  the  processes  are  all  zero  mean  Gaussian,  their  properties  can  be 
completely  specified  by  their  covariance  functions,  i.e., 

E{r(t)r*(u)\H0}  =  Rn(t,u ) 

E{r(t)r*(u)  |.£/i}  =  Rs(t,u)  +  Rn(t,u)  (2) 

where  the  signal  and  clntter  are  assumed  to  be  independent  random  processes  and 

Rs(t,u )  =  E  {s(t)s*(u)} 

Rn(t,u )  =  E  {n(t)n*(u)}  —  N0  S(t  —  u),  (3) 

with  A(-)  denoting  the  Dirac  delta  function.  The  next  step  is  to  expand  the  signal  waveform 
in  terms  of  an  orthonormal  series2: 


s(t)  =  ^sfc0fc(t),  (4) 

k 

where 

sk  =  I  dt  (5) 

and  {4>k(t)}  is  an  orthonormal  set  of  basis  functions  (specihed  later),  i.e., 

/  dt  =  {  J  [”  l  ~  \  (6) 

It  therefore  follows  that  the  signal  waveform  is  completely  specihed  by  its  expansion  coeffi¬ 
cients,  Sk,  which  have  the  properties: 

E{sk}  =  0  (7) 

E{sks^}  =  I J  (f>l(t)Rs(t,u)(j>i(u)  dttdu  (8) 

If  the  orthonormal  basis  functions  also  satisfy 

J  Rs(t,u)<f>i(u )  du  =  A i4>i(t),  (9) 

then 

E  {sfcsz*}  =  A/  5ki-  (10) 

Therefore  are  a  set  of  independent  Gaussian  random  variables  with  zero  mean  and 
variance,  Xk.  The  clutter-only  waveform  can  also  be  expanded  using  the  same  set  of  basis 
functions  so  that 

«(«)  =  E  nk(t>k(t),  (11) 

k 

2For  simplicity,  summations  and  products  written  in  a  form  similar  to  Eq.  (4)  are  indexed  from  k  =  1  to 
oo.  Specific  limits  will  be  used  when  necessary.  Similarly,  integrals  with  no  limits  are  definite  integrals  from 
— oo  to  oo. 
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with 


nk  =  j  <j>*k(t)n(t)  dt 

(12) 

E  {nk}  =  0 

(13) 

E  {nkni  }  =  N0  8ki . 

(14) 

Note  that  the  above  orthonormal  expansion  of  a  random  waveform  whose  basis  functions 
satisfy  the  integral  equation  with  the  covariance  function  as  its  kernel  is  known  as  the 
Karhunen-Loeve  expansion  (KLE)  [3,  4,  5]. 

By  performing  a  Karhunen-Loeve  expansion  on  the  received  signal,  the  binary  detection 
problem  becomes 


H0  ■  rk  =  nk 

Hi  :  rk  =  sk  +  nk  (15) 

where  r(t)  =  J2k  rk(fik(t)  and  rk  =  f  (p*k{t)r{t )  dt.  The  PDFs  of  the  null  and  alternative 
hypotheses  needed  to  calculate  the  likelihood  function  and  hence  the  test  statistic  are  then 
given  by 


f(r(t)\H0)  =  n 

k 

=  n 

k 

tmm  =  n 


l  5 

1  f 

rfc  2 

2nN0  C^P  l 

2  Nq 

1 

-  exp  < 

2vr  (Afc  +  Nq) 

SKP1  _j_ 

2N0  J  v^KVc 


exp 


^{'r  fc}2 

2  Nn 


2  (A k  +  N0)  J  ’ 


(16) 

(17) 


where  K{rfc}  and  Tflrfc}  are  the  real  and  imaginary  parts  respectively.  The  likelihood  function 
is 

L(r(t))  =  Lr(WJ 


No 


n 


(A,  +  N0) 


exp 


E 


A  A 


2N0(Xk  +  No 


■I  rk\‘ 


(18) 


A  detection  is  declared  if  the  likelihood  function  is  greater  than  some  predetermined  decision 
threshold.  It  is  usual  practice  to  use  an  equivalent  but  simpler  form  of  the  decision  function, 
the  test  statistic.  For  this  case,  the  decision  function  is  simply  obtained  by  taking  the 
logarithm  of  the  likelihood  function  and  keeping  only  those  terms  containing  the  received 
signal  so  that  the  test  statistic  is 


T{r(i)}  =  X 

fc 


A  fc 

(A  fc  +  No) 


lrfc|2- 


(19) 
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Since  in  some  practical  situations  where  the  clutter  level  is  unknown,  the  test  statistic  can 
be  further  simplked  by  assuming  the  worst  case  of  when  all  the  signal  components  are  much 
smaller  than  the  associated  clutter  component,  i.e.,  when  <C  Nq.  This  results  in  the  test 
statistic  _ 

7>(i)}  =  £A*M2'  (20) 

k 

which  will  be  adopted  as  the  detector /classifier  structure  for  the  remainder  of  this  report.  A 
schematic  representation  is  shown  in  Fig.  1.  It  can  be  seen  that  this  structure  has  a  simple 
intuitive  interpretation,  i.e.,  the  detection  of  the  target  is  determined  by  a  weighted  sum  of 
the  square  of  the  correlator  outputs,  where  each  target  basis  function  acts  as  a  template  for 
the  correlation.  Note  that  this  structure  is  also  sometimes  referred  to  as  the  separable  kernel 
receiver  [5]  because  of  the  properties  associated  with  the  kernel  of  the  KLE  integral  equation 
(see  Eq.  (9)). 


3  DETECTOR/CLASSIFIER  IMPLEMENTATION 
3.1  Computing  the  Basis  Functions 

The  detector /classifier  structure  derived  in  the  previous  section  is  generic  and  can  be  used 
in  many  different  ways.  This  section  gives  a  detailed  example  of  applying  this  structure  to 
an  echo  location  system  (e.g.,  sonar  or  radar)  as  a  post-processor  for  discriminating  one  type 
of  target  from  the  background,  i.e.,  clutter,  interference  and  other  targets,  etc.  A  summary 
of  the  computational  procedure  required  to  implement  the  post-processor  is  given  as  follows: 

(1)  Generate  the  covariance  matrix  of  the  target  response, 

1  n9 

Rs(i,j )  =  (21) 

9  k=  i 

where  s(ti,  6 *,)  is  the  matched  filter  output  of  the  target  echo,  9f;  is  a  random  variable 
representing  the  aspect  angle  of  the  target  and  Nq  is  the  number  of  aspect  angles 
sampled  3. 

(2)  Solve  the  eigenvalue  problem  (the  discretized  version  of  the  KLE  integral  equation), 

Rs  (hi)  M)  =  ^k4>k{i),  (22) 

3 

where  i,  j,  k  —  1, . . . ,  L.  In  practice,  only  M  of  the  L  most  significant  eigen-solutions 
are  kept.  Of  course,  this  introduces  a  level  of  arbitrariness  since  it  has  to  be  decided 
which  eigen-solutions  are  to  be  discarded. 

3This  assumes  that  the  PDF  of  the  aspect  angles  is  uniform,  which  is  the  most  naive  assumption  that 
can  be  made.  A  more  general  form  of  Eq.  (21)  can  be  made  by  removing  the  factor  1  /Nq  and  multiplying 
each  term  in  the  summation  by  the  probability  of  its  associated  aspect  angle. 
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(3)  Use  the  eigen-solutions,  \k  and  c/>k(t),  to  compute  the  test  statistic 

M 

Tir(t)}  = 

fc= 1 


<t>l(t)r(t)  dt 


(23) 


where  r(t)  represents  the  matched  hlter  output  of  the  received  waveform  for  the  two 
cases:  target  plus  background  and  background  only.  Sample  waveforms  from  these  two 
cases  are  used  to  generate  the  Hq  and  H\  PDFs  of  the  test  statistic.  This  then  allows 
the  decision  threshold  to  be  determined  by  the  Neyman-Pearson  criterion  [3],  i.e., 


T  <  7  Hq 

T>  7  =>  Hi  (24) 


where  the  decision  threshold,  7,  is  set  by  the  desired  false  alarm  probability. 


It  is  evident  from  this  summary  that  information  about  the  target,  the  transmitted  signal 
and  background  are  required  for  the  evaluation  of  the  post-processor’s  performance.  Models 
of  these  three  areas  used  to  demonstrate  the  processing  are  discussed  next. 


3.2  Target  Model 


A  basic  requirement  of  the  above  scheme  is  that  the  target  response  as  a  function  of 
aspect  angle  be  known.  This  information  can  be  obtained  in  several  ways:  held  measure¬ 
ments,  numerical  modelling  or  measurements  using  scale  models.  Each  method  has  its  own 
advantages  and  disadvantages  associated  with  the  computational  and  measurement  difficul¬ 
ties  involved.  For  the  purpose  of  demonstrating  the  technique,  a  simple  mathematical  model 
of  a  spatially  extended  target  is  used.  The  model  consists  of  a  line  of  point  scatterers,  each 
with  a  scattering  amplitude  weighting,  wk,  and  time  delay,  77.  The  important  point  here  is 
that  the  target  exhibits  aspect  angle  dependence.  The  far  held  scattering  response  of  such  a 
target,  using  a  transmitted  waveform  x(t),  is  given  by 


y(t) 


wk(f)  exp  [j2tt f  (t  -  Tk)}  X(f)  df 


(25) 


where 


2  r 

Tk  =  — 
Cm. 


2(k  —  l)d 


cos  9, 


(26) 


r  is  the  range  of  the  target  (assuming  a  monostatic  system),  cm  is  the  speed  of  wave  prop¬ 
agation  in  the  medium,  Ns  is  the  number  of  target  highlights,  d  is  the  spacing  between 
target  highlights,  6  is  the  aspect  angle  of  the  target  measured  relative  to  the  echo-location 
system’s  line  of  sight  to  the  target,  X '(f)  is  the  Fourier  transform  of  x(t),  and  j  =  \f— T.  If 
the  scattering  amplitudes  of  the  individual  highlights  are  independent  of  frequency,  then  the 
target  response  can  be  simplified  to 


Ns 

y(t )  =  ^2wkx(t  -  Tk). 

k=  1 


(27) 


Detection  and  Classification  of  Fluctuating  Targets 


7 


It  then  follows  that  the  matched  filter  output  of  the  received  echo  is 


s(t) 


J  y(£)x*(£  ~  t) 


Ns 


k=  1 


Z({  -  Tk)x*tf  -  t)  (if. 


(28) 


This  then  leads  to  the  requirement  that  the  transmitted  waveform  be  specified  and  hence  a 
discussion  of  the  signal  model. 


To  resolve  the  target  hightlights  and  thus  extract  as  much  information  as  possible  from 
the  target,  the  use  of  waveforms  with  large  time-bandwidth  products  is  desirable4.  A  linear 
frequency  modulation  (LFM)  waveform,  which  is  commonly  used  in  practice,  is  a  particularly 
simple  example  of  a  broadband  signal.  For  mathematical  convenience,  a  LFM  waveform  with 
a  Gaussian  envelope  was  used  in  the  work  reported  here.  In  complex  notation,  the  waveform 
is  given  by 


(  t2  \ 

„  7 tB  9 

C5?  J exp 

j2vr  fct  +  j—t 

1p 

(29) 


where  Tp  is  the  effective  pulse  length  of  the  signal,  fc  is  the  centre  frequency,  and  B  is  the 
bandwidth  of  the  signal.  By  using  Eq.  (29)  in  Eq.  (28)  the  matched  filter  output  of  the 
target  echo  becomes 


Ns 


s(t,  6)  =  TpJ  -  ^  wk  exp  [j2nfc(t  -  Tk )]  exp 


k= 1 


(1  +  (*BT,)»)  2 


2  T2 

p 


(30) 


Note  that  it  is  usually  only  the  complex  envelope  of  the  matched  filter  that  is  of  interest 
because  of  the  uncertainties  associated  with  the  phase  of  the  carrier  signal. 


3.3  Clutter  Model 


In  general,  the  received  signal  will  be  contaminated  by  the  background  environment  so 
that 

r(t,9)  =  s(t,6)  +  n(t),  (31) 


where  n(t)  represents  the  background  held.  In  the  work  reported  here,  It  was  assumed  that 
the  background  is  dominated  by  clutter  whose  behaviour  is  best  modelled  as  stochastic.  In 
particular,  it  was  assumed  that  the  matched  filter  output  of  the  clutter  has  a  phase  angle 
that  is  uniformly  distributed  and  that  its  amplitude  follows  the  Weibull  distribution  [6].  The 
Weibull  PDF  is  given  by 


fz(z) 


cz 

— exp 


(32) 


4There  are  many  waveforms  with  large  time-bandwidth  products  have  approximately  uniform  spectra 
over  a  frequency  interval.  Thus,  the  use  of  such  waveforms  means  that  a  target  can  be  interrogated  with  all 
frequencies  within  the  system  bandwidth.  It  is  in  this  sense  that  such  waveforms  can  be  used  to  extract  the 
maximum  amount  of  information  from  the  target.  It  also  justifies  the  assumption  expressed  in  Eq.  (14). 
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where  z  is  the  value  of  the  random  variable  Z  and  c  and  b  are  the  shape  and  scale  param¬ 
eters  respectively.  The  Weibull  distribution  is  a  convenient  parametrization  of  the  clutter 
amplitude  as  it  can  represent  a  wide  range  of  PDFs  from  Rayleigh  (c  =  2)  to  non- Rayleigh 
by  the  appropriate  choice  of  its  shape  parameter. 


4  PERFORMANCE  EVALUATION 


The  previous  section  described  the  computational  procedure  and  the  various  target,  sig¬ 
nal  and  clutter  models  needed  to  demonstrate  the  workings  of  the  separable  kernel  receiver. 
Its  performance  can  now  be  evaluated  by  choosing  the  following  example  set  of  model  pa¬ 
rameters. 

•  The  target,  composed  of  a  line  of  point  scatterers  (Ns  =  15),  is  divided  into  two  sets: 
(i)  10  scatterers,  each  separated  by  d  —  0.01  m,  with  a  scattering  amplitude  weighting 
Wk  =  1/50  and  (ii)  5  scatterers,  each  separated  by  d  =  10  m,  with  Wk  =  1/5.  The 
propagation  speed,  cm,  is  1500  m/s. 

•  The  transmitted  signal  has  a  carrier  frequency,  fc  —  5  kHz,  effective  pulse  duration, 
Tp  =  0.01  s,  and  bandwidth,  B  =  500  Hz. 

•  The  Weibull  clutter  parameters  are  b  =  l/^/l  +  2/c  (normalised  power)  and  c  ,  where 
c  is  the  free  variable. 


The  target  strength  (the  scattered  energy  in  dB)  and  the  matched  filter  response  are  shown 
in  Figs.  2(a)  and  (b)  respectively.  The  covariance  matrix  of  the  target  used  for  the  determi¬ 
nation  of  the  eigenvalues  and  eigenvectors  of  the  KLE  is  constructed  from  the  matched  filter 
response  of  the  target  at  5  degree  intervals  in  aspect.  The  resultant  normalised  eigenvalues 
are  shown  in  Fig.  2(c)  in  descending  order  of  magnitude  and  the  corresponding  first  five 
eigenvectors  (real  part)  are  shown  in  Fig.  2(d). 

The  performance  of  the  separable  kernel  receiver  can  be  quantified  by  its  ROC  curves, 
specifically,  curves  of  the  probability  of  detection  (Prj)  as  a  function  of  the  probability  of 
false  alarm  (Pfa)  and  signal-to-clutter  ratio  (SCR),  which  is  generically  defined  as 


SCR 


Peak  receiver  response  to  the  target  echo 
Expected  receiver  response  to  the  clutter 


(33) 


Therefore,  PDFs  of  the  test  statistic  were  estimated  by  Monte  Carlo  calculations  because 
analytical  methods  are  generally  unavailable  except  for  certain  special  cases.  Specifically,  a 
modified  “histogram  approach”  was  used  based  on  importance  sampling  [7].  One  thousand 
samples  were  used  for  each  of  the  target-plus-clutter  and  clutter-only  cases.  Furthermore, 
the  target-plus-clutter  PDF  was  obtained  by  sampling  the  target  response  over  all  aspect 
angles.  For  the  studies  presented  in  this  report,  it  was  assumed  that  all  target  aspects  are 
equally  likely,  but  this  need  not  be  the  case  and  the  problem  can  be  easily  modified  to  suit 
other  situations. 
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Examples  of  histograms  of  the  signal-plus-clutter  and  clutter-only  test  statistics  are  shown 
in  Fig.  3  for  a  SCR  equal  to  3  dB  and  9  dB  respectively,  and  c  =  1.75  (non-Rayleigh  clutter). 
A  realisation  of  the  matched  filter  output  for  each  case  is  also  shown.  These  figures  show  that 
the  separation  between  the  signal-plus-clutter  and  clutter-only  PDFs  increase  with  SCR,  as 
expected.  In  particular,  note  that  for  the  case  of  a  SCR  equal  to  3  dB,  the  peak  values  of  the 
clutter  response  are  higher  than  the  peak  values  of  the  target  response.  This  demonstrates 
that  if  detection  is  performed  by  simply  comparing  the  matched  filter  output  to  a  threshold, 
there  is  considerable  risk  that  the  target  would  not  be  detected  and  false  alarms  would  occur. 
Once  the  PDFs  were  found,  the  ROC  curves  were  calculated.  An  example  is  shown  in  Fig.  4 
for  the  target  in  Weibull  clutter  with  c  =  1.75  and  false  alarm  probabilities  Pfa  =  10”2, 
1CT3  and  10-4.  The  idealised  cases  of  Rician  and  Rayleigh  fluctuating  targets  [8]  in  Rayleigh 
clutter  are  also  shown  in  the  same  figure  for  comparison.  The  Rician  case  corresponds 
to  a  non-fluctuating,  dominant  scatterer  in  a  fluctuating  background  of  a  large  number  of 
smaller  scatterers.  The  Rayleigh  case  corresponds  to  a  large  combination  of  scatterers  of 
approximate  equal  strength.  This  case  also  has  the  useful  feature  that  its  ROC  curves  can 
be  calculated  by  mostly  analytical  means.  The  derivation  of  the  PDFs  and  their  detection 
and  false  alarm  probabilities  are  given  in  the  appendix. 

The  differences  between  the  ROC  curves  in  Fig.  4  can  be  explained  by  considering  the 
different  assumptions  used  to  derive  each  detector /classifier.  Matched  filter  receivers  only 
detect  energy  present  in  a  time  series.  Knowledge  of  the  target  and  clutter  characteristics 
is  expressed  only  by  the  PDFs  assumed  or  derived  for  the  target-plus-clutter  and  clutter- 
only  cases.  For  low  SCR  a  non-fluctuating  target  will  almost  never  rise  above  the  clutter, 
and  for  high  SCR  a  fluctuating  target  will  have  some  significant  chance  of  fading  into  it. 
Consequently,  for  low  SCR,  the  ROC  curves  for  the  Rayleigh  fluctuating  target  are  higher 
than  those  for  the  non-fluctuating  (Rician)  target,  but  for  high  SCR,  they  are  below  them.  In 
deriving  the  separable  kernel  receiver,  the  signal  space  of  the  target  is  used,  thus  the  receiver 
looks  for  echoes  that  have  a  specific  “structure.”  Consequently,  it  will  tend  to  ignore  large 
apparent  target  echoes  that  do  not  possess  the  correct  structure,  and  tend  to  detect  echoes 
that  do  even  when  the  SCR  is  low.  Of  course,  the  target  echo  will  fluctuate  because  the 
aspect  is  random.  Therefore,  it  is  reasonable  that  the  separable  kernel  receiver  exhibits  a 
higher  Pd  then  either  matched  filter  receiver  at  low  SCR,  but  does  not  outperform  a  matched 
filter  detecting  a  non-fluctuating  target  at  high  SCR. 

It  is  also  interesting  to  compare  the  performance  of  the  separable  kernel  receiver  with 
that  of  a  peak  detector.  The  peak  detector  uses  no  information  about  the  target  and  simply 
searches  for  the  peak  of  the  matched  filter  output.  Figure  5  shows  the  comparison  between 
the  two  cases,  again  indicating  that  a  proper  use  of  target  information  can  significantly 
improve  the  detection/classification  performance,  which  in  this  case  is  about  a  6  dB  increase 
in  the  SCR  when  Pd  is  equal  to  0.5. 

Next,  consider  how  the  separable  kernel  receiver’s  performance  degrades  when  there  is  an 
adverse  clutter  distribution.  For  a  fixed  Pfa  equal  to  10~3,  Fig.  6  shows  the  ROC  curves  as  a 
function  of  the  Weibull  shape  parameter,  c,  which  controls  the  tail  of  the  clutter  PDF  (c  =  2 
for  Rayleigh,  c  =  1  for  Exponential).  As  in  the  case  for  non- Rayleigh  behaviour,  the  longer 
the  tail  of  the  distribution  (c  decreases),  the  higher  the  chance  of  getting  a  large  amplitude 
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event  that  can  either  mimic  target-like  features  or  mask  the  desired  signal.  Figure  6  shows 
the  detection  performance  is  severely  degraded  as  the  shape  of  the  distribution  becomes 
more  exponential  in  character. 

Until  now  it  has  been  assumed  that  the  target  information  used  to  generate  the  KLE 
eigenvectors  is  fairly  accurate.  In  practice,  this  unlikely  to  be  the  case  because  of  a  number  of 
factors  such  as  measurement  errors,  modelling  and  computational  limitations,  or  distortions 
of  the  target  response  brought  about  by  the  operating  environment.  Figure  7  shows  two 
examples  of  where  the  target  response  (target  strength  and  matched  filter  output)  has  been 
deliberately  altered  from  the  original  response  that  was  used  to  construct  the  KLE  eigenvec¬ 
tors.  The  first  case  consists  of  a  10%  increase  in  the  scatterer  spacing  and  the  removal  of  one 
main  scatterer.  The  second  case  consists  of  a  50%  increase  in  the  scatterer  spacing.  Figures 
8  and  9  show  the  detection  performance  of  these  two  cases  compared  to  the  matched  case. 
As  expected  there  is  a  performance  degradation,  which  is  manifest  as  1.5  to  2  dB  effective 
reduction  in  SCR  for  Pd  =  0.5.  This  implies  that  the  mismatched  separable  kernel  receivers 
will  still  exhibit  a  higher  Pd  than  the  peak  detector.  Stated  alternately,  the  mismatched 
separable  kernel  receiver  will  exhibit  at  least  a  4  dB  effective  increase  in  SCR  over  the  peak 
detector  when  Pd  =  0.5. 


5  DISCUSSION  AND  CONCLUSION 

The  previous  section  (see  Fig.  5)  showed  that  the  separable  kernel  receiver  significantly 
outperformed  the  peak  detector,  i.e.,  a  greater  SCR  is  required  for  the  peak  detector  to 
give  the  same  probability  of  detection  as  the  separable  kernel  receiver.  The  peak  detector 
scheme  in  this  case  is  simply  to  search  for  the  maximum  value  of  the  envelope  of  the  standard 
matched  filter  output.  This  value  is  proportional  to  the  energy  of  the  target  echo.  So  it  is 
pertinent  to  ask  what  is  the  basis  for  this  performance  enhancement  in  the  separable  kernel 
receiver?  To  answer  this,  consider  the  KLE  of  the  signal  and  clutter  waveforms: 


sit)  =  J2skMt) 


n(t)  =  y ^nk(j)k(t). 


The  energy  of  the  signal  waveform  is 


/OO 

\s(t)\2  dt  =  y, 

k= i 


and  similarly  for  the  clutter  waveform,  En  =  Yl'kLi  \nk\2 >  where  the  orthonormal  property  of 
the  KLE  basis  functions  has  been  used.  The  SCR  of  the  peak  detector  is  then  given  by 

rp  sr^oo  |  i2 

SCR0  =  -A  =  ^=l|gfc|2.  (36) 

u  77'  \“^00  II Z  V  / 

^ n  1  ln*M 
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If  the  signal  can  be  sufficiently  characterised  by  the  first  N  terms  of  the  KLE,  i.e., 


S(t)  ~  £  •Sfc^fc  (t)  ! 


then  the  energies  for  the  truncated  versions  of  the  signal  and  clutter  waveforms  are 


Therefore 


ES(N)  =  2~Es  and 

k= 1 
N 

En(N )  =  J]M2. 


SCR(iV)  En  _Er=i  Wt 


SCR0  ~En(N)  Ef=ik|2  ’  V  ' 

so  that  there  is  an  enhancement  in  the  SCR  of  the  truncated  signal  due  to  the  concentration 
of  the  signal  energy  in  a  subset  of  eigenvectors  relative  to  the  clutter. 

Now  consider  the  case  of  the  separable  kernel  receiver  where  its  outputs  for  the  signal 
and  clutter  components  are 


Ttm  =  £  Afc|sfc|2  and 

k=  1 
M 

T  {n(t)}  =  ^  A k  Wk\2  ■ 


Then  define  the  SCR  of  the  separable  kernel  receiver  output  as 


SCRr  = 


__  T  {»(«)} 
"  T  {«(*)} 


so  that 


SCRr  Ehi  kl2  E2Li  Kl! 


QPP  v^oo  I  |2  \  |  1 2  *  V  7 

Let  the  eigenvalues,  A*,,  be  normalised  to  their  largest  value  and  placed  in  descending  order. 
If  the  signal  energy  is  essentially  confined  to  a  finite  number  of  eigenvectors,  as  indicated  by 
the  decay  of  the  eigenvalues,  then 


£a 


k  I  '-’fc  —  /  I &k 


and  by  the  same  argument  as  above, 


(46) 
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In  effect,  the  eigenvalues  provides  a  softer,  less  abrupt  truncation  of  the  signal  expansion 
than  would  otherwise  be  the  case  for  a  simple  termination  of  the  series.  However,  the 
enhancement  in  the  SCR  in  both  cases  comes  from  the  compact  nature  of  the  signal  energy 
when  the  signal  is  decomposed  according  to  the  KLE  basis  functions. 

In  conclusion,  this  report  presents  a  technique  for  the  detection  and  classification  of 
complex  targets  whose  echoes  are  strongly  sensitive  to  the  target’s  aspect  dependence.  It  has 
been  shown  that  a  priori  information  about  the  target  can  be  used  to  enhance  the  detection 
and  classification  performance.  The  KLE  was  used  to  form  the  basis  of  a  receiver  structure 
that  exploited  the  compactness  of  the  signal’s  energy  relative  to  that  of  the  background. 
In  addition,  the  KLE  has  the  advantage  of  providing  a  “natural”  signal  space  basis  into 
which  the  received  waveform  can  be  decomposed.  This  signal  decomposition  was  achieved 
by  solving  the  eigenvalue  problem  associated  with  the  KLE  of  the  target  signal.  The  resultant 
eigenvectors,  weighted  by  their  eigenvalues,  played  the  role  of  templates  for  the  detection 
and  classification  of  the  target. 
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APPENDIX 


PROBABILITY  DENSITY  FUNCTION 
OF  THE 

SEPARABLE  KERNEL  RECEIVER  OUTPUT 


This  appendix  outlines  the  derivation  of  the  probability  density  functions  (PDFs)  and 
their  corresponding  detection  and  false  alarm  probabilities  of  the  separable  kernel  receiver 
for  the  case  of  a  fluctuating  Rayleigh  target  in  Rayleigh  clutter.  Recall  that  the  test  statistic 
for  the  separable  kernel  receiver  is 


M 

Tir(t)}  =  ^2^k\rk\2 

k= 1 


(Al) 


and  so  the  aim  is  to  obtain  the  PDFs  of  T  for  both  the  target-plus-clutter  and  clutter  only 
cases.  To  begin,  assume  that  the  in-phase  ^{ry})  and  the  quadrature  (S{ry})  components 
of  the  received  waveform  follow  a  zero  mean  Gaussian  distribution  with  the  variance  given 
by 

H0  :  a l  —  N0  (clutter  only) 

Ih 


=  A*,  +  N0  (target  plus  clutter). 


(A2) 


Then  the  PDF  of  the  amplitude  of  ry  is 


tn  i  \  I  Del  \  \rk\ 2 

/(W)  =  ^ex  p 


(A3) 


It  is  convenient  to  introduce  a  change  of  variables  so  that  the  test  statistic  is  T  =  J2kLi  uk-> 
where  uk  =  Xk  \rk\~ ■  The  PDF  of  the  new  random  variable  then  becomes 


9  M  = 


exp 


2A  ka2k 

and  its  corresponding  characteristic  function  for  uk  is 


uk  \ 


(A4) 


*uk(t)  = 


g(uk)e3Ukt  duk 


1 


1  -  j2A  ko\t' 


(A5) 
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where  j  =  \J  —  1.  Therefore,  the  characteristic  function  for  T  is 


M 


n  *-.(*) 


E 


Ak 

l-  j2\k<j2kt' 


(A6) 


where  Ak  are  the  residues  of  the  partial  fraction  expansion  of  the  product  (assuming  no 
degeneracy  in  the  poles).  By  taking  the  inverse  Fourier  transform  of  the  characteristic 
function,  the  PDF  of  the  test  statistic  is 


M 


!t(T)  =  Y. 

k=  1 


Ak 

2A^feXP 


T 


(A7) 


and  the  false  alarm  and  detection  probabilities  are  given  by 


Pt(T) 


fr(t)  dt 


IT 

M 


1V1  s 

Y  Ak  exp 

k=\  '■ 


T 

2A  kal 


(A8) 


for  each  of  the  two  cases  given  in  Eq.  (A2).  This  last  set  of  equations  allows  the  detection 
threshold  to  be  determined  for  a  given  false  alarm  probability.  This  may  be  accomplished 
using  Newton’s  method.  The  threshold  can  also  be  used  to  calculate  the  probability  of 
detection. 


Eigenvalue  (normalised) 
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Test  Statistic ,  T, 
compared  to 
Neyman-Pearson 
Threshold 


y[KK(t) 

Figure  1:  The  structure  of  the  separable  kernel  receiver. 


(a)  Target  Strength  (dB) 
90  40 


Eigenvalue  Number 


(b)  Matched  Filter  Output  (dB) 


Aspect  angle  (deg) 
(d) 


Figure  2:  Some  properties  of  the  target  response:  (a)  target  strength  as  a  function  of  aspect, 
(b)  matched  filter  output  as  a  function  of  aspect,  (c)  the  eigenvalues  of  the  KLE,  and  (d) 
the  corresponding  first  five  eigenvectors  of  the  KLE. 


Matched  Filter  Amplitude 
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(a) 

Signal  and  Clutter  Distributions:  peak  SCR  =  3  dB 


(b) 

Signal  and  Clutter  Distributions:  peak  SCR  =  9  dB 


Example  of  a  MF  output:  peak  SCR  =  9  dB 


Figure  3:  The  histograms  of  the  signal-plus-clutter  and  clutter-only  test  statistics  for  two 
cases  of  SCR:  (a)  3  dB  and  (b)  9  dB.  An  example  of  a  matched  filter  output  realization  is 
also  shown  for  each  case. 
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Figure  4:  Three  sets  of  ROC  curves  show  the  performance  of  the  separable  kernel  receiver 
for:  (solid  lines)  the  target  in  Weibull  clutter  (c  =  1.75),  (broken  lines)  a  Rician  fluctuating 
target  in  Rayleigh  clutter,  and  (dash-dot  lines)  a  Rayleigh  fluctuating  target  in  Rayleigh 
clutter.  Within  each  set,  each  curve  corresponds  to  a  false  alarm  probability  of  10~2,  10-3 
and  1CT4  from  left  to  right. 
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c  =  1.75 


Figure  5:  The  ROC  curves  of  the  separable  kernel  receiver  (solid  lines)  and  the  peak  detector 
(broken  lines)  for  falsm  alarm  probabilities  of  ICO2,  ICO3  and  10-4  in  Weibull  clutter  (c  = 
1.75). 


^=io~3 


Figure  6:  The  ROC  curve  of  the  separable  kernel  receiver  as  a  function  of  the  Weibull  shape 
parameter  for  a  false  alarm  probability  of  ICO3. 
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(a) 


(b) 


Figure  7:  The  target  strength  and  matched  filter  output  as  a  function  of  aspect  for  two  cases: 

(a)  the  spacing  between  the  main  scatterers  is  increased  by  10%  with  one  less  scatterer;  and 

(b)  the  spacing  between  scatterers  is  increased  by  50%. 


Figure  8:  The  ROC  curves  of  the  separable  kernel  receiver  showing  the  degradation  in 
performance  when  there  is  a  10%  mismatch  in  the  scatterer  spacing  and  one  less  main 
scatterer. 
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10%  mismatch  in  scatterer  spacing  and  one  less  scatterer 


Figure  9:  The  ROC  curves  of  the  separable  kernel  receiver  showing  the  degradation  in 
performance  when  there  is  a  50%  mismatch  in  the  scatterer  spacing. 


